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ABSTRACT 

A  two-dimensional  shear-deformable  beam  element  is  developed  in  this  investigation  for 
the  analysis  of  the  large  rotation  and  large  deformation.  Using  the  absolute  nodal 
coordinate  formulation  and  a  continuum  mechanics  approach,  the  assumption  of  Euler- 
Bemoulli  and  Timoshenko  beam  theories  are  relaxed.  The  effect  of  the  shear  deformation 
is  accounted  for  without  the  need  for  introducing  Timoshenko’s  shear  coefficient.  By 
using  the  absolute  coordinates,  the  nonlinear  strain-displacement  relationships  are  used  to 
define  a  relatively  simple  expression  for  the  elastic  forces,  while  the  mass  matrix  of  the 
beam  remains  constant.  As  a  consequence,  the  centrifugal  and  Coriolis  forces  are 
identically  equal  to  zero.  Surprisingly,  the  more  general  model  developed  in  this 
investigation  leads  to  a  significant  saving  in  computer  time  as  compared  to  non-shear 
deformable  models  presented  in  previous  investigations. 
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1.  INTRODUCTION 

In  Euler-Bemoulli  beam  theory  it  is  assumed  that  the  cross-section  remains  rigid  and 
perpendicular  to  the  neutral  axis  of  the  beam  [1, 13].  In  this  theory,  the  effect  of  the  shear 
deformation  is  neglected.  In  Timoshenko  beam  theory  on  the  other  hand,  the  cross- 
section  does  not  remain  perpendicular  to  the  beam  neutral  axis.  Nonetheless,  the  cross- 
section  remains  rigid.  A  shear  coefficient  is  introduced  in  order  to  account  for  the  shear 
deformation  [12].  In  this  investigation,  a  two-dimensional  shear-deformable  beam 
element  based  on  the  non-incremental  absolute  nodal  coordinate  formulation  [10]  is 
developed.  In  this  approach,  only  absolute  coordinates  and  global  slopes  are  used  to 
define  the  element  nodal  coordinates  without  the  need  for  using  infinitesimal  or  finite 
rotations.  With  the  definition  of  the  nodal  coordinates  in  the  global  coordinate  system,  no 
transformation  is  required  to  determine  the  element  inertia  or  elastic  forces.  Using  this 
coordinate  representation  with  the  appropriate  element  shape  function,  exact  modeling  of 
the  rigid  body  dynamics  can  be  achieved.  Using  the  non-incremental  absolute  nodal 
coordinate  formulation,  the  resulting  mass  matrix  of  the  finite  element  is  a  constant 
matrix  and  the  centrifugal  and  Coriolis  forces  are  identically  equal  to  zero.  Another 
advantage  of  using  the  absolute  nodal  coordinate  formulation  is  its  simplicity  in 
formulating  the  generalized  forces  and  imposing  some  joint  constraints. 

A  problem  encountered  in  the  implementation  of  the  non-incremental  absolute 
nodal  coordinate  formulation  is  the  formulation  of  the  elastic  forces.  Shabana,  et  al  [6,  8, 
9,  10,  11]  proposed  two  methods  for  formulating  the  elastic  forces  of  the  two- 
dimensional  beam  element.  In  the  first  method,  a  local  element  coordinate  system  is 
introduced  for  the  convenience  of  describing  the  element  deformation.  This  approach 
leads  to  a  complex  expression  for  the  elastic  forces  even  when  a  linear  elastic  model  is 
used.  In  the  second  method  [3]  a  continuum  mechanics  approach  is  used  to  obtain  the 
elastic  forces  without  introducing  the  local  element  coordinate  system.  In  this  continuum 
mechanics  approach,  nonlinear  strain-displacement  relationships  must  be  used  in  order  to 
obtain  zero  strain  under  an  arbitrary  rigid  body  motion.  Nonetheless,  the  previous  models 
developed  using  the  continuum  mechanics  approach  are  based  on  Euler-Bemoulli  beam 
theory  which  dose  not  account  for  the  shear  deformation.  It  is  the  objective  of  this 
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investigation  to  develop  a  model  for  the  elastic  forces  for  two-dimensional  beam 
elements  that  accounts  for  shear  deformation  by  using  a  general  continuum  mechanics 
approach  without  introducing  a  local  element  coordinate  system.  This  new  model  relaxes 
the  assumptions  of  Euler-Bemoulli  and  Timoshenko  beam  theories  and  does  not  require 
the  use  of  a  shear  coefficient. 

This  paper  is  organized  as  follows.  In  section  2,  the  displacement  field  and  the 
kinematic  equations  of  the  new  two-dimensional  shear-deformable  beam  element  are 
presented.  The  constant  mass  matrix  of  the  beam  element  that  accounts  for  the  rotary 
inertia  is  defined  in  section  3.  In  section  4,  the  formulation  of  the  elastic  forces  using  the 
general  continuum  mechanics  approach  is  discussed.  The  formulations  of  the  generalized 
external  forces  and  moments  associated  with  the  element  nodal  coordinates  using  the 
virtual  work  are  presented  in  section  5.  In  section  6,  the  equations  of  motion  are 
presented  and  the  use  of  the  proposed  model  in  the  analysis  of  the  large  deformation  of 
simple  and  multibody  mechanical  systems  is  demonstrated.  Section  7  presents  a  summery 
and  the  conclusions  drawn  from  the  study  presented  in  this  paper. 

2.  KINEMATICS  OF  THE  SHEAR-DEFORMABLE  BEAM 

For  the  two-dimensional  shear-deformable  beam  element,  the  displacement  field  is 
defined  in  the  global  coordinate  system  as 


a^+a^x+ajy+a^xy+a^x2  +a5x3 
b0 +blx+b2y+b3xy+b/ix1  +b5x3 


(1) 


where  r,  as  shown  in  Fig.  1,  is  the  global  position  vector  of  an  arbitrary  point  on  the  beam 
element  cross-section,  at  and  b,  are  the  polynomial  coefficients,  and  x  and  y  are  the  spatial 
coordinates  defined  in  a  beam  coordinate  system.  The  spatial  coordinate  x  is  chosen  to  be 
along  the  beam  axis  (0  <*</),  where  l  is  the  element  length.  Note  that,  the  assumed 

displacement  field  depends  on  y  in  order  to  account  for  the  shear  deformation. 
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In  the  simplified  Euler-Bemoulli  beam  theory,  the  effect  of  the  shear  deformation 
is  neglected  [1].  The  basic  assumption  in  the  simplified  Euler-Bemoulli  beam  theory  is 
that  the  cross-section  of  the  beam  remains  normal  to  the  beam  neutral  axis  as  shown  in 
Fig.  2a.  The  beam  cross-section  at  any  point  along  the  beam  neutral  axis  can  be  defined 
by  the  Frenet  frame.  The  Frenet  frame  has  one  of  its  axis  tangent  to  the  beam  neutral  axis 
and  the  other  axis  perpendicular  to  the  beam  neutral  axis  [7].  The  tangent  vector  t_  can  be 
defined  by  dr/dx .  In  a  shear-deformable  beam  model,  the  cross-section  of  the  beam  does 
not  remain  normal  to  the  neutral  axis,  as  shown  in  Fig.  2b.  As  a  result,  the  tangent  to  the 
neutral  axis  cannot  be  used  to  define  the  cross-section.  In  order  to  demonstrate  that  the 
shape  function  of  Eq.  1  accounts  for  the  shear  effect,  consider  an  arbitrary  vector  Ar, 
which  is  defined  in  the  beam  cross-section  as  shown  in  Fig.  3.  Using  the  displacement 
field  defined  by  Eq.  1,  it  can  be  shown  that 

Ar  =  r-r  0  =y~~  (2) 

dy 

where  r  is  the  global  position  vector  of  an  arbitrary  point  P  in  the  cross-section  with 
coordinates  (x,y),  and  rv=0  is  the  position  vector  of  the  corresponding  point  P0  on  the 

beam  centerline  with  coordinates  (x,0).  The  preceding  equations  show  that  any  arbitrary 
vector  drawn  on  the  beam  cross-section  can  be  defined  by  the  vector  dr  /dy,  and  as  a 
consequence  the  vector  dr/dy  defines  the  cross-section  of  the  beam.  In  the  absolute 
nodal  coordinate  formulation,  the  global  position  vector  of  an  arbitrary  point  on  the  beam 
can  be  written  as 


r  = 


r\ 


=  Se 


ri 


(3) 


where  S  is  the  global  element  shape  function,  and  e  is  the  vector  of  nodal  coordinates. 
The  vector  of  the  element  nodal  coordinates  e  is  given  by 

e  =  h  <?2  e3  e4  e5  e6  ei  %  e\o  en  en]T  (4) 
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The  vector  of  nodal  coordinates  includes  the  global  displacements 


e\  ~  rl|jc=0  5  e2  ~  r2  |x-0  5  e7  ~r\  x=/  3  eS  ~  r2  I  x=l 


and  the  global  slopes  of  the  element  nodes  that  are  defined  as 
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The  element  shape  function  S  must  have  a  complete  set  of  rigid  body  modes  that  describe 
arbitrary  rigid  body  translational  and  rotational  displacements.  Using  the  element  nodal 
coordinates  given  by  Eq.  4,  the  element  shape  function  can  be  defined  as  follows 

^  0  Is  2  0  ls3  0  s4  0  ls5  0  ls6  0 

o  —  ( 0 ) 

0^0  ls2  0  ls3  0  $4  0  ls5  0  ls6  J 

where  the  functions  s,  =s(^,rj)  are  defined  as 

*=  l-3^2  +  2^3,  s2=S-2ijl+Si , 

s4=3f2-2£3,  s5=-42+l?, 

and  ^  =  y ,  77  =  y ,  /  is  the  element  length. 

3.  CONSTANT  MASS  MATRIX 

The  global  position  vector  of  an  arbitrary  point  on  the  shear-deformable  beam  is  given  by 
Eq.  3.  By  differentiating  this  equation  with  respect  to  time,  the  absolute  velocity  vector 
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can  be  defined  as  r  =  Se.  This  vector  can  be  used  to  define  the  kinetic  energy  of  the 
element  as 


T  =  -  [prTrdV  =  — e7 
2i  2 


1  • 
e  =  — e  M„e 


(6) 


where  V  is  the  volume,  p  is  the  mass  density  of  the  beam  material,  and  Ma  is  the  mass 
matrix  of  the  element.  The  mass  matrix  in  Eq.  6  is  constant  and  symmetric.  Using  the 
shape  function  given  by  Eq.  5,  the  mass  matrix  is  given  by 


M.  =  JpS1 
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(7) 


where  m  is  the  total  mass  of  the  finite  element,  /  is  the  element  length,  J\  is  the  first 
moment  of  mass  defined  as  Jx=^pydV ,  and  J2  is  the  second  moment  of  mass  defined  as 


J2  =  \py2dV. 
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4.  ELASTIC  FORCES 


In  this  section,  the  nonlinear  strain-displacement  relations  are  used  to  develop  an 
expression  for  the  elastic  forces  of  the  beam  element.  The  deformation  gradient  can  be 
defined  as  [2] 


dr, 

dr\_ 

dx 

dy 

dx 

dr2 

dr2 

dx 

dy. 

S2xe  S2ye 


(8) 


where  Sa  =  — ,  Siv  =  — and  S ,•  is  the  i  th  row  of  the  element  shape  function.  The 

dx  y  dy 

Lagrangian  strain  tensor  sm  can  be  written  as 


8 


m 


(jrj-i)=i 

"erSae-l 

erSce 

V  /  2 

.  erSce 

erSAe-l 

(9) 


where  I  is  the  identity  matrix,  Sa  =S1rcSu+S2.tS2;c,  SA  =  s[yS,_l,  +S2>,S2>, ,  and 
Sc  =  s[xSly  +  S2xS2y .  It  should  be  noted  that  strain  tensor  is  symmetric  thus  it  can  written  in 
a  vector  form  as 

e  =  [e,  s2  e3]r  (10) 

where  £|  =  ^-(erSae-l)  ,  e2  =  ^-(erSAe-l) ,  and  e3  =  i-erSce  . 

A  general  expression  for  the  strain  energy  can  be  written  using  the  strain  vector  e 
and  the  stress  vector  a  =  [a1  <j2  ffj]7  as  follows  [4,  5] 


Using  the  constitutive  equations,  the  stress  vector  is  related  to  the  strain  vector  by 


cj  =  Ee  (12) 

where  E  is  the  matrix  of  the  elastic  constants  of  the  material.  For  isotropic  homogenous 
material,  the  matrix  E  can  be  expressed  in  terms  of  Lame’s  constants  X  and  fj.  as 


X+2fi 


X  0 
X+2(i  0 

0  2  n 


(13) 


where  X= - ,  u  = - ,  E  is  the  Young’s  modulus  of  elasticity,  and  vis  the 

(l  +  v)(l-2v)  2(1 +  v)  &  J 

Poisson’s  ratio  of  the  beam  material.  Using  Eqs.  12  and  13,  the  strain  energy  can  be 
rewritten  as 


U  =  -  JV'EsrfF 

V 


The  vector  of  the  elastic  forces  Qe  can  be  defined  using  the  strain  energy  U  as 


(14) 


de 


(15) 


where  K  is  the  stiffness  matrix  which  can  be  written  as 


K  =  ( A+2yu)  Kj  +  A  K2  +  2//  K3 


(16) 


where 


K,  =  ^  J[Sfll  (erSae  -  1)+SM  (erSAe- l)]dV , 

V 

^i=\  j[sa,  (erSAe  -  1)+SM (erSae - 1  )]dV , 
H  v 
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K3  =^-  jsd(erSae)  dV , 

*  V 

Sai  =  +  STa ,  Sbl  =  Sb  +  si ,  and  Scl  =  Sc  +  . 

Note  that  the  general  expression  for  the  elastic  forces  obtained  in  this  investigation  using 
the  continuum  mechanics  approach  and  the  nonlinear  strain-displacement  relationships  is 
much  simpler  than  the  expression  obtained  in  previous  investigations  [6,  9]  using  the 
element  local  coordinate  system  and  linear  strain-displacement  relationships.  The  general 
expression  obtained  in  this  paper  automatically  captures  the  effect  of  geometric 
centrifugal  stiffening. 

5.  FORMULATION  OF  THE  GENERALIZED  EXTERNAL  FORCES 

The  virtual  work  can  be  used  to  develop  the  vector  of  the  generalized  external  forces 
[10].  The  virtual  work  due  to  an  externally  applied  force  F  acting  on  an  arbitrary  point  on 
the  element  is  given  by 


SWe  =  Fr<5r  =  Fr  S  de  =  Q^<Se 


(17) 


where  r  is  the  position  vector  of  the  point  of  application  of  the  force,  and  Qe  =  S7  F  is  the 

vector  of  the  generalized  forces  associated  with  the  element  nodal  coordinates.  For 
example,  the  virtual  work  due  to  the  distributed  gravity  of  the  finite  element  can  be 
obtained  using  the  shape  function  and  Eq.  17  as 


J[0  -pg]S<5e  =  -mg 
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12 


0  0  0  - 
2 


12 


0  0 


<5e 


(18) 


This  defines  the  vector  of  generalized  distributed  gravity  forces  as 


Q  =  -mg 


0-0  —  000-0  00 
2  12  2  12 


-I  T 


(19) 
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When  an  external  moment  M  acts  on  the  cross-section  of  the  beam,  the  virtual 
work  due  to  this  moment  is  given  by 


SWM  =MSy  (20) 

i 

where  y  is  the  angle  of  rotation  of  the  cross-section.  The  orientation  of  a  coordinate 
system  attached  to  the  cross-section  can  be  defined  using  the  following  transformation 
matrix 


cos(y) 
sin  (y) 


-sin(y) 

cos(y) 


dr2 
1  dy 
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-dt\ 

dy 


£ r± 

dy 
dr2 
d y 


/  =  . 
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dr2 

\dyj 


(21) 


which  yields 


cos(y)= 


1  dr, 
f dy ’ 


.  /  X  1  drx 

sin  y  )= - - 

K)  fdy 


(22) 


Using  these  two  equations,  it  can  be  shown  that  the  virtual  change  in  the  cross-section 
orientation  angle  can  be  defined  as 


<5y= 


di\ 

_^L  s 

fdj£ 

dy 

Kdyj 

f2 


(23) 


If  a  concentrated  moment  is  applied  for  example  at  node  A  as  shown  in  Fig.  4,  the 
generalized  force  vector  due  to  this  moment  is  defined  as, 


Q  M  - 


oooo 


Me,  -Me 


fl 


n 


^-  0  0  0  0  0  0 


(24) 


9 


6.  NUMERICAL  EXAMPLES 


Using  the  obtained  expressions  of  the  elastic  and  inertia  forces,  the  equations  of  motion 
can  be  obtained  using  the  absolute  nodal  coordinates  as 

Mae  =  Q  (25) 

where  Q  is  the  vector  of  generalized  nodal  forces  including  the  elastic  force  vector  Qe 
and  external  force  vector  Qi ,  and  e  is  the  vector  absolute  accelerations.  As  previously 
mentioned,  the  centrifugal  and  Coriolis  force  vectors  are  equal  to  zero  since  the  mass 
matrix  is  constant. 

In  this  section,  two  examples  are  considered  in  order  to  demonstrate  the 
performance  of  the  proposed  beam  model.  The  two  examples  considered  are  the  free 
falling  of  a  flexible  pendulum  under  its  own  weight,  and  a  four  bar  mechanism  with  a 
flexible  connecting  rod. 

The  first  example  considered  is  the  free  falling  two-dimensional  beam  shown  in 
Fig.  5.  The  beam  is  connected  to  the  ground  by  a  pin  joint  at  one  end.  The  beam  has 
length  of  1.2  m,  circular  cross-sectional  area  of  0.0016  m2,  second  moment  area  8.533E- 
06,  a  mass  density  of  5540  kg/m3,  Poisson’s  ratio  of  0.3,  and  a  modulus  of  elasticity  of 
0.700E  06.  The  beam  is  initially  horizontal  with  zero  initial  velocity  and  free  to  fall  under 
the  effect  of  gravity.  Two  cases  are  considered  in  this  example.  In  the  first  case,  the 
gravity  constant  is  assumed  to  be  9.81  m/s2  while  in  the  second  case,  the  gravitational 
acceleration  is  increased  to  50  m/s  .  The  simulations  of  the  beam  are  performed  using 
different  number  of  elements.  Figure  6  shows  the  position  of  the  tip  point  of  the  beam 
using  6  and  12  finite  elements  when  the  gravity  constant  is  equal  to  9.81  m/s2.  It  is  clear 
from  the  results  presented  in  this  figure  that  there  is  a  good  agreement  between  the  two 
models.  The  results  demonstrate  that  the  solution  converges  with  small  number  of 
elements. 

Since  the  falling  pendulum  is  conservative  system,  the  sum  of  the  energy 
components  must  remains  constant,  that  is 
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(26) 


'Yri+Ui+Vi=const 

i 

where  t‘  is  the  element  kinetic  energy,  U  is  the  element  strain  energy,  v'  is  the  element 
potential  energy,  and  ne  is  the  number  of  elements  of  the  system.  Figure  7  shows  the 
energy  balance  for  6-element  model.  The  results  presented  in  this  figure  shows  that  the 
total  energy  remains  constant.  The  results  obtained  using  the  6-element  model  are  the 
same  as  the  results  of  the  12-element  model.  Figure  8  shows  the  displacement  of  the  tip 
point  of  both  models  under  gravitational  acceleration  of  50  m/s2.  Very  small  differences 
between  the  two  models  can  be  observed.  The  energy  balance  obtained  using  the  6- 
element  model  under  the  gravity  constant  of  50  m/s2  is  shown  in  Fig  9.  The  same  results 
are  obtained  using  the  12-element  model.  Figure  10  shows  the  large  deformation 
configurations  of  the  falling  beam  at  different  time  steps  under  gravitational  acceleration 
of  50  m/s2  using  48  finite  elements. 

The  shear-deformable  model  developed  in  this  investigation  relaxes  the 
assumption  of  Euler-Bernoulli  beam  theory.  In  this  model,  the  cross-section  does  not 
remain  perpendicular  to  the  beam  centerline  due  to  the  effect  of  the  shear  deformation. 
Figure  1 1  shows  a  comparison  between  the  results  obtained  using  this  new  model  and  the 
Euler-Bernoulli  beam  model  previously  presented  by  Berzeri  et  al  [3]  who  used  a  shape 
function  which  does  not  depend  on  y.  The  results  obtained  in  this  figure  are  obtained 
using  12  elements.  It  is  important  to  point  out  that  the  shear-deformable  model,  because 
of  the  dependence  of  the  shape  function  on  y,  leads  to  a  simpler  expression  for  the  elastic 
forces.  More  surprisingly,  significant  saving  in  computer  time  was  achieved  using  the 
shear-deformable  model.  It  was  observed  that  the  shear-deformable  model  is  two  times 
faster  than  the  Euler-Bernoulli  model.  The  results  presented  in  Fig.  1 1  shows  that  there  is 
a  very  good  agreement  between  the  shear-deformable  beam  model  and  the  model  based 
on  the  Euler-Bernoulli  beam  theory  since  a  thin  beam  is  used.  When  the  cross-section 
changed  to  0.001  m2  and  second  moment  area  1.593E-06  m4,  noticeable  differences  are 
observed  as  shown  in  Fig.  12. 

The  shear-deformable  model  developed  in  this  investigation  also  relaxes  the 
assumption  of  the  Timoshenko  beam  theory  since  it  allows  for  the  plane  deformation  of 
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the  cross-section,  that  is  the  cross-section  remains  plane  but  not  rigid.  The  deformation  of 
the  cross-section  can  be  measured  by  the  deviation  of  the  norm  of  the  vector  dr  /  dy  from 
one.  Figure  13  shows  the  norm  of  this  vector  at  node  5  and  node  9  as  a  function  of  time  in 
the  case  of  free  falling  pendulum  and  gravity  constant  of  50  m/s2. 

The  second  example  considered  in  this  numerical  study  is  the  four  bar 
mechanism.  The  mechanism  consists  of  a  rigid  crankshaft,  a  flexible  connecting  rod,  and 
a  flexible  follower.  The  initial  configuration  of  the  mechanism  is  shown  in  Fig.  14.  The 
geometric  and  inertia  properties  of  the  four  bar  mechanism  are  shown  in  Table  1.  The 
crankshaft  is  driven  by  a  moment  defined  by  the  function  shown  in  Fig.  15.  In  this 
simulation,  one  element  is  used  for  the  crankshaft,  6  elements  for  the  connecting  rod,  and 
4  elements  for  the  follower.  The  longitudinal  and  transverse  deflection  of  the  mid  point  of 
the  connecting  rod  is  shown  in  Fig.  16.  Also,  for  this  multibody  example  the  solution 
converges  with  small  number  of  elements. 

7.  SUMMARY  AND  CONCLUSIONS 

In  this  investigation,  a  shear-deformable  beam  model  based  on  the  non-incremental 
absolute  nodal  coordinate  formulation  is  developed  for  the  large  rotation  and  large 
deformation  analysis.  The  beam  model  leads  to  exact  modeling  of  the  rigid  body 
dynamics  and  leads  to  zero  strain  under  an  arbitrary  rigid  body  displacement. 
Furthermore,  the  model  relaxes  the  assumptions  of  Euler-Bemoulli  and  Timoshenko 
beam  theories.  The  cross-section  in  the  new  model  does  not  remain  rigid  and  does  not 
remain  perpendicular  to  the  beam  centerline.  By  using  a  continuum  mechanics  approach, 
the  model  leads  to  a  relatively  simple  expression  for  the  elastic  force. 

While  the  model  accounts  for  the  effect  of  the  rotary  inertia  and  shear 
deformation,  the  mass  matrix  remains  constant.  As  a  consequence,  the  centrifugal  and 
Coriolis  forces  are  identically  equal  to  zero.  Two  numerical  examples,  a  free  falling 
pendulum  and  a  multibody  four  bar  mechanism,  were  used  to  demonstrate  the  use  of  the 
new  beam  model.  Numerical  results  obtained  in  this  study  demonstrated  that  the  solutions 
obtained  using  the  new  model  converges  much  faster  as  compared  to  the  non-shear- 
deformable  models  presented  in  previous  investigations.  Furthermore,  significant  saving 
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in  computer  time  is  achieved  by  using  the  more  general  shear-deformable  model  that  is 
based  on  the  nonlinear  strain-displacement  relationships. 
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Table  1.  The  geometric  and  inertia  properties  of  the  four  bar  mechanism 


Body 

m  (kg) 

A  (m2) 

/(m4) 

/(m) 

E  (Pa) 

V 

Crankshaft 

0.681 

1.257E-03 

1.257  E-07 

0.200 

1.000  E+09 

0.3 

Coupler 

2.474 

1.960  E-03 

3.068  E-07 

0.900 

5.000  E+06 

0.3 

Follower 

1.470 

7.068  E-0 

3.976  E-08 

0.519 

5.000  E+08 

0.3 
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a)  The  beam  in  the  undeformed  configuration 


Fig.  1.  The  position  vector  of  an  arbitrary  point  on  the  beam  cross-section 


a)  The  Euler-Bemoulli  beam  without  shear  deformation 


dr 


b)  The  beam  element  with  shear  deformation 


Fig.  2.  The  beam  deformation  assumptions 
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Fig.  3.  The  position  vector  of  an  arbitrary  point  on  the  beam  cross-section 


Fig.  4.  Externally  applied  moment 


Vertical  Displacement  (m) 


Fig.  9.  Energy  balance  of  the  beam  ( g  =  50  m/s2) 


X  (m) 


Fig.  10.  The  falling  flexible  pendulum  at  different  time  steps  using  48  elements 

( g  =  50  m/s2 ) 
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Fig.  11.  Comparison  between  the  results  Euler-Bemoulli  beam  model  and  shear 

deformable  beam 


Fig.  12.  Comparison  between  the  results  Euler-Bemoulli  beam  model  and  shear- 
deformable  beam  for  thick  cross-section 


Fig.  15.  The  driving  moment  applied  to  the  crankshaft 


Fig.  16.  The  longitudinal  and  lateral  deformation  of  the  mid  point  of  the  coupler  link 


